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In this ab initio study, we focused on the influence of chemical composition and correlation effects 
on the electronic and magnetic properties of Sr 2 FeMoC >6 (SFMO). We considered the appearance of 
the two most important lattice defects - oxygen vacancies and antisite defects. Different DFT-based 
methods were used, combining the advantages of a multiple scattering Green’s function method, a 
spin-polarized fully relativistic linear-muffin-tin-orbital method, and the Vienna ab initio simulation 
package. Correlation effects were included with a self-interaction correction and a Hubbard U 
method. Using the magnetic force theorem, we determined the magnetic exchange coupling strength 
in SFMO with and without defects and determined the Curie temperatures Tc via Monte Carlo 
simulations. We found a ferrimagnetic ground state and verified its stability against electronic 
correlations or defects. The total magnetic moments decreased for both kinds of defects, either 
vacancies or antisite defects. This agreed well with the experimental observations. On the other 
hand, we could predict a strong variation of Tq as a function of the defect concentration. Tc could 
be increased by adding more oxygen vacancies while it could be reduced with increase of antisite 
disorder. The x-ray absorption spectra and x-ray magnetic circular dichroism at the Fe and Mo 
Ij 2 , 3 , Fe K and Mo M 2,3 edges were simulated and compared with the available experimental results. 

The spectra agreed only after considering also oxygen vacancies. Thus, experimental observations 
of mixed valence states of the Fe ions in SFMO might be attributed to a certain amount of oxygen 
vacancies. 

PACS numbers: 61.72.-y 75.30.Et 75.47.Lx 78.20.-e 
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I. INTRODUCTION 

Double perovskites A 2 BB'Oe (A = alkaline earth or 
rare earth and BB' are heterovalent transition metals 
such as B = Fe, Cr, Mn, Co, Ni; B' = Mo, Re, W) 
often demonstrate intrinsically complex magnetic struc¬ 
tures and a wide variety of physical properties (see [ I ] for 
a review article on these materials). 

Several experimental and theoretical studies have 
demonstrated that the double perovskite system 
Sr 2 FeMoC >6 (SFMO) and other related materials exhibit 
a ferrimagnetic (FiM) half-metallic ground state with a 
high Curie temperature of 420 K [2]. The physical origin 
of the magnetoresistance in SFMO is half-metallicity [3], 
i.e., the material is an insulator in one of the spin chan¬ 
nels, but a metal in the other. This leads to a complete 
spin polarization at the Fermi level, which immediately 
suggests their application as a source of spin polarized 
charge carriers in spintronic devices. Therein, SFMO 
will be used mainly as a thin film, and many attempts 
were made to grown high quality SFMO films [4-8]. But 
all these films yield a reduced Curie temperature, which 
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is up to 80 K smaller than for bulk samples [7, 8], and 
even for bulk SFMO the theoretical magnetic moment 
of 4 /xb was never experimentally observed. Such varia¬ 
tions were attributed to lattice defects like grain bound¬ 
aries or point defects. The latter includes swapping of Fe 
and Mo ions, antisite disorder (ASD), and oxygen vacan¬ 
cies (Vo)- Several studies contributed to a theoretical 
understanding of these point defects and their effect on 
the magnetic properties [9-15] but only few of them ad¬ 
dressed explicitly the magnetic exchange interactions and 
the Curie temperature, Tq [9, 10, 14]. Thus, we studied 
from ab initio the variation of the magnetic properties of 
SFMO with point defects deploying a Green’s function 
(GF) method with the coherent potential approximation 
(CPA) in order to understand experimental differences in 
the magnetic moments and the Curie temperatures. We 
calculated magnetic coupling constants using the mag¬ 
netic force theorem [16]. They are used in a second step 
to calculate the Curie temperature. The whole proce¬ 
dure provides usually correct magnetic properties of sys¬ 
tems with itinerant and strongly localized electrons, if 
the crystalline structure is known and the applied den¬ 
sity functional approximation describes adequately the 
electronic structure of the studied system [17—23] . 

Before including defects, we compared our numerical 
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results of the electronic structure and magnetic proper¬ 
ties of defect-free SFMO with previous theoretical [9- 
15] and experimental studies as a first validation step. 
Experiments on SFMO included numerous methods like 
photoemission spectroscopy (PES) [24-26], Mossbauer 
spectroscopy [ 7, 28], neutron scattering [29-31] or x-ray 
absorption spectroscopy (XAS) and x-ray magnetic cir¬ 
cular dichroism (XMCD) measurements [32-35]. Since 
the XAS and XMCD spectroscopy is very sensitive to 
the electronic structure and the local environment, fit¬ 
ting calculated spectra to the available experimental re¬ 
sults can provide important information about the chem¬ 
ical composition, the ionic valency and the degree of elec¬ 
tronic correlations in the system. Thus, we simulated the 
XAS and XMCD spectra with our method and found a 
better agreement than previous density functional calcu¬ 
lations [36]. Our results indicate as well a mixed valency 
of the Fe ion (Fe 2+ or Fe 3+ ) observed by experiments 
[28, 29, 33, 34], but we can also follow from our study 
that the mixed valency might be mainly caused by oxy¬ 
gen vacancies. 

In contrast to some previous theoretical studies [10, 
24], we obtained half-metallic SFMO only if electron 
correlation effects were considered - in accordance with 
[11, 15, 37]. Within the adapted experimental lattice 
structure [38], the expected FiM state of SFMO was the 
stable magnetic ground state as long as the electronic 
structure showed half-metallicity. We traced the grad¬ 
ual transition from strong band magnetism in metallic 
SFMO to the ferrimagnetism in half-metallic SFMO and 
the corresponding critical temperatures by varying the 
electronic correlation parameter U e g. Finally, we mod¬ 
eled point defects up to 15 % and 10 % for ASD and oxy¬ 
gen vacancies. The calculated magnetic moments and 
Tq showed significant variations in good agreement with 
available comparable studies but were also strongly de¬ 
pending on the correct description of the underlying elec¬ 
tronic structure. 

In the following section, we describe the lattice struc¬ 
ture of SFMO used in our calculations. Details about 
the calculation techniques are given in the next section. 
The influence of the electronic correlations and the lat¬ 
tice defects on the electronic structure is discussed in 
Sec. IV. Their effect on the magnetic properties, such as 
the magnetic moment or the Curie temperature, is shown 
in Sec. V. We present in the Sec. VI the XAS and XMCD 
calculations for the SFMO compound and compare them 
with previously measured spectra. The results are sum¬ 
marized in the last section. 


II. LATTICE STRUCTURE 

For our calculations of SFMO, we adopted the exper¬ 
imentally found double perovskite structure, where the 
oxygen atoms provide an octahedral environment around 
the Fe and Mo sites. The FeOg and MoOg octahedra al¬ 
ternate along the three cubic axes, while the Sr atoms 


TABLE I. Variation of the experimental lattice constants and 
oxygen positions of SFMO from several references. Fe oc¬ 
cupies consistently (0, 0, 0). The Wyckoff positions of the 
oxygen atoms are (0,0, z) and (x,y, 0). The positions of the 
other atoms are described in the text. 


a (A) 

c (A) 

c/a/y/2 

X 

y 

z 

Ref. 

5.5729 

7.9077 

1.003353 

0.248 

= X 

0.235 

[40, 44] 

5.573 

7.902 

1.002612 

0.2464 

= X 

= X 

[38] 

5.5752 

7.89251 

1.001009 

0.2610 

0.2333 

0.2474 

[41] 

5.587 

7.894 

0.999087 

0.248 

= X 

0.235 

[45] 


occupy the hollow site formed by the corners of the FeC>6 
and MoO e octahedra at the body-centered positions (see 
Fig. 1). 

SFMO was found to be cubic (Fm3m) in the param¬ 
agnetic phase, but changes into a tetragonal-type struc¬ 
ture below a critical temperature [31, 39]. A representa¬ 
tive choice of experimentally observed lattice constants is 
collected in Tab. I. Their variation between the smallest 
and largest value is quite small (0.25% and 0.17% for 
a and c, respectively) and the cubic symmetry would be 
represented by a c/a ratio of \/2. Thus, we see from the 
data in Tab. I that the deviation from the cubic symme¬ 
try is small as well. 

In general, the Sr atoms occupy the 4c? Wyckoff po¬ 
sitions (0, i/ 2 , 1 / 4 ). Depending on the publication, the Fe 
atoms occupy the 2 a Wyckoff pos. (0, 0, 0) and Mo the 4c? 
Wyckoff pos. (0,0, Y 2 ), or vice versa. The positions of the 
oxygen atoms is not definite and varies between different 
studies (see Tab. I). They sit either exactly between Mo 
and Fe, which gives the body centered tetragonal struc¬ 
ture type (/4/mmm; No. 139 [38, 40]) or the oxygen 
octahedra are slightly distorted (J4/m, No. 87 [41, 42]). 
In any case, there are two types of oxygen atoms with O- 
at the 4e Wyckoff pos. (0, 0, z) and O xy at the 8 h Wyck¬ 
off pos. (x, y, 0) while the Fe atom position is (0,0,0) 
(see Fig. 1). 

Since the changes of the lattice constants are small, 
we used for simplicity the more symmetric body centered 
tetragonal structure type and the lattice constant from 
[38] as input for our study. As long it is not stated other¬ 
wise, we used the primitive unit cell with one functional 
unit of SFMO (see blue arrows in Fig. 1). 

III. COMPUTATIONAL DETAILS 

For the microscopic understanding of the SFMO com¬ 
pound, we combined the theoretical results of three 
computational methods, namely the multiple scattering 
Green’s function method HUTSEPOT [46], the spin- 
polarized fully relativistic linear-muffin-tin-orbital (SPR- 
LMTO) method [47-49] and the Vienna ab initio simu¬ 
lation package (VASP) [50, 51]. The main investigation 
of the electronic and magnetic properties were conducted 
with HUTSEPOT, whereas the XAS and XMCD spectra 
were calculated within the SPR-LMTO, and necessary 
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FIG. 1. (Color online) The double perovskite structure of 
SFMO. The colored polyhedra visualize the octahedral sur¬ 
roundings of the Fe and Mo atoms (orange and blue) [13]. 
Following from the tetragonal symmetry, two different oxygen 
positions appear (marked with O xy and O z ). The tetragonal 
supercell is shown by the black solid lines and the red ar¬ 
rows. It contains two functional units with two Fe sites (Fel, 
Fe2). The black dashed lines and the blue arrows indicate the 
primitive unit cell. 


structure relaxations were achieved with VASP. 

a. Electronic structure The orbital resolved den¬ 
sity of states (DOS) were calculated with HUTSEPOT 
within the full-charge density approximation to the crys¬ 
tal potential. For the Green’s function, the angular mo¬ 
mentum cutoff was Z max = 3 and the Brillouin zone inte¬ 
gration was done with a fc-point mesh of 12 x 12 x 12. The 
complex energy contour was integrated over 24 Gaussian 
quadrature points. With the Green’s function G(E ) of 
the system, all quantities of interest like the local density 
of states (LDOS) follow in a straightforward way. 

The exchange-correlation functional of a GGA-type 
was used very successful in the context of SFMO [3, 52]. 
Within the Green’s function method, we used in partic¬ 
ular the version of Perdew, Burke and Ernzerhof (PBE) 
[53, 54]. In order to deal with the strong correlations 
in the oxide system, the self-interaction correction (SIC) 
[ 16] and the GGA+t/ approach of Dudarev et al. [55] 
implemented in the Green’s function method were used 
( U e g = U — J). For the SIC, only the LDA exchange- 
correlation functional can be used. 

We note, that the calculation of the electronic struc¬ 
ture in SFMO is very sensitive to the choice of the atomic 
radii and the crystalline structure. We used the atomic 
radii, obtained from the largest overlap of atomic poten¬ 
tials. Using smaller radii for the Fe atoms led to a half- 
metallic solution independent of the exchange-correlation 
functional - either LDA or GGA. We observed the same 
also within the SPR-LMTO method. 


b. Magnetic properties We calculated the inter¬ 
atomic exchange coupling parameters Jij from ab ini¬ 
tio using the magnetic force theorem [16] implemented 
within HUTSEPOT. The underlying magnetic reference 
state for these calculations was the FiM solution - Fe 
moments align parallel to each other and Mo moments 
align antiparallel to the Fe moments. This was also iden¬ 
tified as stable magnetic ground state by the calculation 
of the corresponding magnon spectrum. On the other 
hand, the magnetic moments of Mo are often considered 
as induced moments. Such moments will vanish above 
the Tq and alter, thereby, the numerical results for Tq. 
We can model a paramagnetic (PM) state at high tem¬ 
perature via the disordered local moment (DLM) theory 
[56, 57] and the corresponding results will be compared 
with the conventional calculations. Often, this procedure 
allows a better estimation of Tq [23]. 

In any case, the J can be used in a Monte 
Carlo (MC) simulation to obtain the critical tempera¬ 
tures (details in [18, 58, 59]). Therefore, the unit cell in 
Fig. 1 was repeated 20 x 20 x 20 times to form a large 
cluster. We considered also periodic boundary condi¬ 
tions and restricted the calculation only to the magnetic 
atoms of the unit cell. The starting point was a high- 
temperature disordered state above the critical temper¬ 
ature Tq. In the course of the simulations, the temper¬ 
ature was stepwise reduced until magnetic ordering was 
reached. For each temperature T, the thermal equilib¬ 
rium was assumed to be reached after 20 000 MC steps. 
The thermal averages were determined over 20 000 addi¬ 
tional MC steps. Tq was then obtained from the tem¬ 
perature dependency of the magnetic susceptibility and 
the heat capacity within a numerical uncertainty range 
of ±5K. For comparison, Tq was also calculated with 
the mean-field approximation (MFA) and random-phase 
approximation (RPA) [60]. 

c. Simulation of x-ray spectra The x-ray absorp¬ 
tion and dichroism spectra were calculated with the SPR- 
LMTO method taking into account the exchange split¬ 
ting of the core levels. The approach for the XAS and 
XMCD simulations is described in our previous papers 
[C ]. The finite lifetime of a core hole was accounted for 
by folding the spectra with a Lorentzian. The widths of 
core level spectra A. 2 = 1.14eV and /l 3 = 0.41 eV for 
Fe, T h2 = 1.83 eV, T u = 1.69 eV and T M23 = 2.1 eV for 
Mo were taken from [62]. The finite apparative resolu¬ 
tion of the spectrometer was accounted for by a Gaussian 
of width 0.6 eV. 

In our SPR-LMTO simulations, the basis set con¬ 
sisted of the s, p, and d, LMTO’s for the Sr, Fe, Mo 
and O sites. The fc-space integrations were performed 
with the improved tetrahedron method [63] and the self- 
consistent charge density was obtained with 1063 irre¬ 
ducible fc-points. We used a relativistic generalization 
of the LDA+t/ method, which takes into account the 
spin-orbit coupling so that the occupation matrix of lo¬ 
calized electrons becomes non-diagonal in spin indexes. 
This method is described in detail in our previous pa- 
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per [64] including the procedure to calculate the screened 
Coulomb U and exchange J integrals, as well as the Slater 
integrals F 2 , E 4 , and F e . The constrained LSDA calcu¬ 
lations produce J = 0.85 eV for the Fe site in SFMO. 
In order to be consistent with the other methods, the 
Hubbard U has been used as an external parameter. 

d. Defects We investigated the influence of two 
types of point defects in Sr 2 [Fei_ a: Mo 3 ,][Moi_ a; Fe x ] 06 - 5 , 
namely, oxygen-deficiency 6 and the ASD x. The latter 
means a swap of an Fe ion with a Mo ion and vice versa 
keeping the stoichiometry of SFMO constant. 

Both were simulated with the coherent potential ap¬ 
proximation in the GF method, while the lattice struc¬ 
ture has to be kept static. Although some experimental 
studies report no random distribution but segregation 
into small clusters and formation of antiphase boundaries 
[65], the ASD was modeled by interchanging randomly Fe 
and Mo. At the low concentration limit, the CPA is ap¬ 
propriate to reproduce the electronic structure of both 
randomly distributed and clustered disorder. Advanced 
schemes [66], which go beyond the single-site approxima¬ 
tion of the CPA need still high computational resources 
for larger unit cells and are beyond the scope of this inves¬ 
tigation. For the oxygen vacancies, a certain percentage 
of empty spheres was introduced at the lattice sites of the 
oxygen ions. The typical oxygen-deficiency 5 ranges be¬ 
tween 0.006 to 0.36 [67-69]. This represents 0.1 at. % to 
6 at. % of the total oxygen amount in defect-free SFMO. 

In addition, we used a supercell approximation with 
a cell of two functional units of SFMO (see red arrows 
in Fig. 1). This supercell was relaxed with VASP includ¬ 
ing either ASD or a single oxygen vacancy close to Fe2. 
For the antisite defects, we found, similarly as in [11], no 
significant relaxation. The relaxation with VASP for a 
single oxygen vacancy led to a distance of 1.9722 A be¬ 
tween Fe2 and the vacancy - an increase by 1.5%. The 
correspondent distance between the Fel ion and the oxy¬ 
gen vacancy was reduced by 0.5% to 4.4054 A. 


IV. ELECTRONIC STRUCTURE 

Although the half-metallic solution was previously ob¬ 
tained within the LDA and the GGA with different cal¬ 
culation techniques [3, 70], the Green’s function method 
yielded a metallic solution similar as in [11, 15, 37] with 
a nonzero density of states in both spin channels at the 
Fermi energy, Ep, independent of the use of LSDA or 
GGA (see Figs. S2 and S3 [71]). We also crosschecked 
the calculation of the DOS with the other methods and 
found similar results. 

From this metallic state, a band gap in the major¬ 
ity spin channel emerges when correlation corrections 
are applied (either GGA+t/ [11, 37] or SIC [72]), and 
a half-metallic state forms [71]. Optical measurements 
determined a band gap of 0.5 eV, which should appear 
between the Mo t\ & states and the Fe e| states [73]. This 
gap is rather small considering the position of the Fe 


e| states (« 1.3 eV) measured by photoemission spec¬ 
troscopy (PES) [24-26]. It might be that intermediate 
energy levels were introduced by defects, and allow addi¬ 
tional optical transitions. Hence, we expected the Fe ej 
[74] states roughly at 1.3 eV below Ep. Applying the SIC 
in our calculations, the localization of the Fe d states 
was strongly overestimated and the states were too far 
below Ep , although the half-metallic ground state is well 
represented within the SIC [72]. In order to obtain a bet¬ 
ter agreement for the Fe e| state with the experimental 
observations, we calculated the DOS for several correla¬ 
tion parameters Ujg applied only to Fe d states [71], and 
found a good agreement for Ujg = 2 eV (see dashed line 
in Fig. 2). The main changes occur for the Fe states at 
Ep , which crossed Ep for t/JJf 1 eV and the band gap 



FIG. 2. (Color online) The upper panel shows the experi¬ 
mental and simulated PES of SFMO taken as a copy of Fig. 
5(a) of Ref. [24] (some labels were removed). The lower 
panel represents the DOS for SFMO calculated with the GF 
method and = 2eV. The gray area shows the total DOS 
for defect-free SFMO. The colored regions indicate the LDOS 
for the d states of Fe (reddish) and Mo (bluish). The green 
dotted DOS includes 8 at. % of randomly distributed oxygen 
vacancies. Both plots are scaled with respect to the same en¬ 
ergy scale and the Fermi energy at zero (black vertical line). 
The dashed line indicates the position of the Fe e| in the 
experimental PES. 
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emerged [71]. Then, almost all eft (eft) states of Fe are 
occupied (unoccupied), which is represented by a formal 
Fe 3+ valence state. For the Mo ions, almost all states 
are unoccupied, except the delocalized t\ states, which 
could be considered as a Mo 5+ state. This situation was 
found in all self-consistent solutions with GGA+/7 and 
represents the generally accepted ground state of SFMO 
[3, 11, 14, 72]. Thus, we use it as well in the description 
of the magnetic properties below. 

Although the finding of an appropriate value for Ujg, 
which was described above, is a rather standard pro¬ 
cedure when dealing with strongly correlated systems, 
several aspects might influence the comparison between 
theoretical DOS and experimental PES data. Here, we 
emphasize two aspects. On the one hand, the calculated 
Fe e| states had a bandwidth of roughly 1 eV (see Fig. 2), 
whereby the PES measurement can not define the exact 
position of the Fe e| states. On the other hand, the con¬ 
clusion about the contribution of the states to the PES 
is based on numerical calculations of defect-free SFMO 
[24]. Saitoh et al. [24] noticed possible ASD in their sam¬ 
ples (up to 10 at. %) but could not account for defects in 
their theoretical interpretation. When we include within 
the CPA some atomic percent of ASD inside SFMO, the 
DOS is strongly varied as observed by many other stud¬ 
ies [10-12, 15] . In particular, some states of the major¬ 
ity spin channel appear also at Ep and reduce the spin 
polarization, which is also observed experimentally [75]. 
In contrast, a low concentration of oxygen vacancies does 
not alter the half-metallic character of SFMO [ 1, 12, 76], 
but additional Fe (ft and Mo dft states become occupied 
[11]. This indicates a lower valency of the Fe and Mo ions. 
Electrons of the removed oxygen atoms occupy states of 
Fe and Mo ions. A similar qualitative behavior was ob¬ 
served within our calculations of oxygen-deficient SFMO 
with CPA (green dashed line in Fig. 2). The intensity of 
the states in the spin down channel right below the Fermi 
energy is increased as obtained in Ref. [11], and the Fe 
e| states are shifted to lower energies. In addition, the 
averaging within the CPA broadens all states (see Fig. 2). 

Taken as a whole, a direct comparison between the 
DOS and the PES data is complicated when considering 
all the various influences of electron correlations or pos¬ 
sible defects. In order to have a point of reference for 
later discussions, we believe that Ujg = 2eV is a suitable 
compromise considering all available data. It is applied 
in the following, if not stated otherwise. 


V. MAGNETIC PROPERTIES 

The observed transition from a metallic to a half- 
metallic state influences as well the magnetic properties. 
We investigated theoretically the magnetic moments, the 
magnetic exchange coupling, and the Curie temperature 
of SFMO in dependence on the electron correlations and 
defect concentrations. All three numerical methods men- 



FIG. 3. Magnon spectrum of half-metallic SFMO (Uftg = 
2eV). Two well separated modes can be identified. One of 
them (positive) is associated mainly with the Fe magnetic 
atoms, while the other one (negative) with the Mo moments. 


tioned above showed a good qualitative agreement de¬ 
spite of the different exchange functionals. The ferri- 
magnetic ground state was always energetically favor¬ 
able when comparing possible magnetic configurations 
and we did not observe a hint to spin spiral ordering 
as reported earlier [10]. In order to verify the magnetic 
ground state, we calculated the magnon spectrum for 
half-metallic SFMO (see Fig. 3). The spectrum shows 
two well separated modes. In accordance with our anal¬ 
ysis of the corresponding site-resolved magnonic Green 
function, the positive mode is mainly associated with 
the collective precession of Fe magnetic moments, while 
the negative mode reflects the properties of Mo magnetic 
moments. The different signs of the modes indicate on a 
ferrimagnetic ground state, which is confirmed by our to¬ 
tal energy calculations. The different observations might 
follow from our choice of taking the experimental lattice 
structure of SFMO instead of the undisturbed cubic lat¬ 
tice structure [10]. However, Solovyev [10] found that the 
FiM ground might be stabilized including breathing dis¬ 
tortions of the oxygen octahedra. If the oxygen ion came 
closer to the Fe ions, the proposed spin spiral ordering 
became less favorable with respect to the FiM state but 
the half-metallicity was gone. These observations match 
ours. Our calculations within the GGA yielded for the 
experimental structure (where dFe-o/^Fe-Mo < 0.5, see 
Tab. I) a metallic ground state as well. 


A. Magnetic moments 

We calculated the spin and orbital moments in the 
FiM ground state. The Fe spin and orbital moments 
are parallel, whereas the spin and orbital Mo moments 
are antiparallel to each other, in accordance with Hund’s 
third rule (see Tab. II). Our calculations for the metallic 
solution (weak electron correlation regime) provide for Fe 
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a lower spin magnetic moment than in earlier calculations 
(3.72 hb to 3.97/is [3, 36, 77-79]). As soon as we observe 
half-metallic SFMO (Ujg > leV), the Fe spin moment 
of ss 4 is in a good agreement with the experimental 
neutron diffraction measurements [30] and the theoretical 
ideal value for the Fe 3+ /Mo 5+ or Fe 2+ /Mo 6+ valency 
configuration (see Tab. II). It remains quite stable and 
insensitive to the strength of the electron correlations. 
The orbital Fe magnetic moment of 0.041 /kb is very close 
to the FPLMTO result of Jeng and Guo [79] (0.043/re) 
and the LMTO results of Kanchana et al. [36] (0.042 /kb). 
The calculated spin and orbital moments of the Mo atom 
are always larger in comparison with earlier calculations 
[3, 36, 77-79]. 

We note that the spin and orbital magnetic moments 
were obtained from the XMCD experiments by using the 
sum rules, which relate the integrated signals over the 
spin-orbit split core edges of the circular dichroism to the 
ground state orbital and spin magnetic moments [S ]. It 
is well known that the application of the sum rules some¬ 
times results in an error up to 50% [81]. To investigate 
the possible error of the sum rules in the case of SFMO we 
compare the spin and orbital moments obtained from the 
theoretically calculated XAS and XMCD spectra through 
sum rules with directly calculated LDA+t/ values. The 
number of the transition metal 3d electrons is calculated 
by integrating the occupied d local density of states inside 
the corresponding atomic spheres, which gives the values 
n Fe = 6.737 and n^° = 4.552. The sum rules reproduce 
the spin magnetic moments within 4 % and 8 % and the 
orbital moments within 21% and 14% for Fe and Mo, 
respectively (see Tab. II). 

The lower reported spin moments of the XMCD ex- 


TABLE II. The calculated spin m s and orbital m 0 magnetic 
moments (in /kb) of Fe and Mo in defect-free SFMO compared 
with the measured magnetic moments. The moments m B and 
m 0 were obtained with the SPR-LMTO (1). In addition, we 
show m s determined with the GF method (2). The sum rules 
are applied for the theoretically calculated XMCD spectra 
in the LDA+17 approximation (SPR-LMTO). Experimental 
uncertainties are given in brackets behind the values. 


method 

m s 

Fe 

m 0 

m s 

Mo 

m 0 

LSDA (1) 

3.503 

0.046 

-0.456 

0.034 

LDA+t/ (1) 

4.104 

0.041 

-0.504 

0.047 

Sum rules^ 1 - 1 

4.283 

0.052 

-0.464 

0.041 

GGA< 2 ' 

3.090 


-0.555 


GGA+[/ (2) 

3.921 


-0.525 


Ref. [ 3] a 

3.05(20) 

0 .02(2) 

-0.32(5) 

-0.05(5) 

Ref. [35] a 

2.80(30) 

0.093(10) 

-0.36(3) 

-0.037(15) 

Ref. [30] b 

4.1(1) 


0 .0(1) 


Ref. [36] c 

3.72 

0.042 

-0.29 

0.020 

Ref. [77] c 

3.97 


-0.39 



a XMCD measurement 
b Neutron diffraction experiment 
c Theory 


oxygen content 6 — 5 
5.4 5.55 5.7 5.85 6 



ASD (at. %) oxygen content (at. %) 


FIG. 4. Total magnetic moment m to t in dependence on 
(a) the ASD, or (b) oxygen deficiency 5 obtained with the 
GF method for metallic, nearly half-metallic and half-metallic 
SFMO (GGA, = 0.75 eV, and = 2eV, respectively). 
Results of previous references are included. 

periments compared to the theoretical value for Fe 3+ 
(4 /kb) were always attributed to antisite disorder [9, 33, 
67, 82]. We used again the CPA to investigate randomly 
distributed defects. The introduced Fe ions at the Mo 
sites (FeMo) coupled strongly antiferromagnetically to the 
intrinsic Fe ions (see Fig. 7). Indeed, the total mag¬ 
netic moment of SFMO is linearly reduced for ASD up 
to 15at.%, independent of the strength of the electron 
correlations (see Fig. 4(a)). The slope of this reduction 
is in good agreement with the estimations by previous 
Monte Carlo simulations [9] and experiments [82]. 

On the other hand, oxygen vacancies reduce also 
the total magnetic moment [67], which agrees well with 
our theoretical results within the CPA. The reduction is 
purely linear for half-metallic SFMO {U^g = 2eV) (see 
Fig. 4(b)). The offset for 6 = 0 can be attributed to 
ASD of roughly 5 at. % to 8 at. % [ 7]. We want to em¬ 
phasize the change of the slope for Ufg = 0.75 eV, where 
we observe the same slope as that of Ujg = 2eV only 
after adding few atomic percent of oxygen vacancies. As 
observed for the DOS, the randomly distributed oxygen 
vacancies have a similar effect as the correlation param¬ 
eter Ufg. Both reduce the electron hopping between the 
Fe and Mo orbitals, and considering the oxygen vacancies 
results in half-metallic SFMO with the correct magnetic 
moment (compare blue and red lines in Fig. 4(b)). 


B. Magnetic exchange interactions and Curie 
temperatures 

In contrast to the magnetic moment, the calculated 
magnetic exchange parameters between the magnetic 
atoms and the critical temperatures are more sensitive 
to changes in the DOS. The most prominent coupling 
constants of the order of several meV have only a very 
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FIG. 5. (a) Sketch of the orientation of the important mag¬ 
netic exchange interactions Joi' X between Fe and Mo ions 
(i is the shell index, X = Fe or Mo) [ ]. (b) The Jot~ X 

in the FiM or PM (DLM) reference state in dependence on 
Ul fj. (c) Curie temperatures calculated with MFA, RPA, or 
MC method. They are compared with the experimental Tc 
range. The observed ground state change in the DLM model 
is marked with open and full circles (see text). The result of 
the calculation with SIC is shown for comparison. 

restricted range up to 7.9 A. This includes only the in¬ 
teractions up to the next nearest neighbor Fe ions J F |~ Fo 
(sketched in Fig. 5(a)). As it is discussed below, only 
these magnetic couplings dominate the magnetic behav¬ 
ior with respect to the electron correlations. The more 
distant exchange constants (up to 12.49 A) were one or¬ 
der of magnitude smaller and became even more reduced 
with increasing Ujg. All the others could be considered to 
be zero. Due to the tetragonal structure, all magnetic ex¬ 
change interactions show a small asymmetry with respect 
to those with a component in 2 direction. For clarity, only 
the coupling constants in the cc-y-plane are shown. 


We considered two different calculation schemes: 

(i) On the one hand side, the calculations of were 
performed in the FiM reference state (red and blue in 
Fig. 5(b)). We observe a strong antiferromagnetic (AFM) 
coupling between Fe and Mo ions. Its absolute size be¬ 
comes reduced around Ujg « leV and increases again 
afterwards. The nearest neighbor Fe-Fe coupling con¬ 
stants, J F i“ Fe , decrease linearly from 4meV to ImeV 
with increasing U^g. It might be correlated to the lin¬ 
ear changes in the Fe e| states. This behavior is typi¬ 
cal for the application of U e g. The stronger localization 
of the orbitals leads to a decrease of the electron hop¬ 
ping and, thereby, to a decrease in the magnetic coupling 
strength. Nevertheless, the behavior of the next near¬ 
est Fe-Fe coupling constants, J F |“ Fe , is more complicated 
(see Fig. 5(b)). These show a strong linear reduction un¬ 
til the coupling switches to a negative value (AFM). At 
Ujg « 1 eV, the slope is reversed but the J F |~ Fo coupling 
constants remain negative. Only the absolute values de¬ 
crease. 

These results deviate from the former calculations 
of magnetic exchange interactions of Solovyev [10]. He 
found for purely cubic SFMO no metallic DOS and, 
hence, no positive J F |' Fe but a rather strong AFM cou¬ 
pling of — 21.5eV. The observations might qualitatively 
match our results for half-metallic SFMO - positive 
J F f" Fe and negative J F 2 ~ Fc - but together with a small 
Fe-Mo coupling of —1.21eV, the FiM state is not stable 
according to his criterion of instability [10] 

2 (j 0 F r Fe + JS? Fe ) - 4 >r Mo < 0. (i) 

The FiM stability could also not be improved by shift¬ 
ing the unoccupied Fe states up by A U and thereby 
increasing the intra-atomic exchange splitting [10]. We 
assume that the lack of stability of the FiM state is caused 
by too large J F | _Fe and too small J F f" Mo constants previ¬ 
ously calculated. When we apply the criterion (1) to our 
magnetic exchange constants, the left hand side is always 
larger than zero, and the FiM state is stable. Neverthe¬ 
less, the qualitative trend observed for increasing A U [1( )] 
can be confirmed by our calculations including Ufg, if we 
compare only half-metallic SFMO (Ufg > 1 eV) - reduc¬ 
tion of the strength of </ F f“ Fc and J F f~ Fe and increasing 
strength of J F f" Mo . 

(ii) In a second set of calculations, we estimated 
the exchange coupling within the DLM approach, which 
models the PM state high above the Curie temperature. 
In this case, the induced moments vanish. The obtained 
magnetic coupling constants (now all J Fe ~ Mo vanish) are 
also shown in Fig. 5(b). They have a similar tendency 
as the magnetic exchange interactions obtained from the 
FiM reference state, which were discussed at first. 

These were considered to calculate the Curie 
temperature, using either the mean-field approximation 
(MFA), the random-phase approximation (RPA), or the 
Monte Carlo (MC) method. For the latter, the cou¬ 
pling constants up to a distance of 12.49 A were used. 
All three methods show, in particular for half-metallic 













SFMO, qualitatively a similar non-linear variation with 
increasing U F § (see Fig. 5(c)). Since the MFA usually 
overestimates Tq and RPA and MC method have similar 
results, we limit our discussion in the following to the 
results of the MC simulations. 

The Jij obtained in the FiM reference state within 
the GGA yielded a Tq of ss 604K, i.e., 200 K above the 
measured values. With increasing Ujg but still a metal¬ 
lic DOS, the Curie temperature decreases linearly un¬ 
til Ujg = 1 eV. This corresponds with the sign change 
of Jq 2 ~ Fb and the decrease of For half-metallic 

SFMO, the theoretical critical temperatures increased 
slightly but stayed always 200 K below the experimen¬ 
tal Tc (compare shaded area in Fig. 5(c)). The increase 
is observed despite the decrease of Joi’ Fe - The stronger 
AFM coupling between Fe and Mo sites mediates an ad¬ 
ditional FM coupling and increases Tq. We obtain also 
from the MC simulations a FiM ground state. For the 
sake of completeness, the T^ c for a SIC calculation is 
plotted as a green point and matches with the Tc in the 
limit of large U F g (see Fig. 5). 

We note the importance of the magnetic moments at 
the Mo site, when considering the DLM model. The van¬ 
ished Mo moments led to an increase of the magnitude 
of the nearest neighbor Fe-Fe coupling constants. Above 
Ujff = 1 eV, the first and second neighbor interactions be¬ 
come equal in magnitude but with different sign. Such a 
behavior alters the critical temperature dramatically and 
we observed only low ordering temperature and a ground 
state with vanishing averaged saturation magnetization 
(non-magnetic (NM), full circles in Fig. 5(c)). Although 
the DLM is expected to give a better description of the 
magnetic phase transition, its possible failure is not com¬ 
pletely new. The induced moments of Mo might remain 
also above the transition temperature and the current 
single-site approximation of the local moments might not 
cover those correctly. An example of this issue was found 
during the investigation of the magnetic properties of Ni 
where only a non-local DLM approach could describe the 
moments sufficiently [83, 84). 

In summary, the theoretical description underesti¬ 
mates Tc of half-metallic SFMO with respect to the ex¬ 
perimental results. The value of Ujg = 2 eV appears as a 
good compromise between half-metallic DOS or high Tc- 
However, there might be various reasons for the reduced 
value - defects like ASD or oxygen vacancies, additional 
electron doping or shortcomings in the description of the 
electronic structure. Their effect was estimated in the 
following sections. 


C. The Curie temperature and defects 

Using the CPA, we simulated at first the effect of ASD 
and oxygen vacancies up to 15 at. % or 10at.%, respec¬ 
tively. Previous simulations or experimental measure¬ 
ments are only available for ASD [9, 10, 14, 86, 87]. They 
saw often a slight reduction of Tc with ASD [9, 85, 8 j 


oxygen content 6 — 5 



FIG. 6. Curie temperatures in dependence on the (a) ASD 
and (b) oxygen deficiency 5 obtained with the MC method for 
metallic, nearly half-metallic and half-metallic SFMO (GGA, 
Ugff = 0.75 eV, and U^g = 2eV, respectively). Results of 
previous references are included. The FiM order is lost for 
more than 5 at. % ASD (marked by dashed line). 
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FIG. 7. Magnetic exchange interactions Jij for (a) metallic 
and (b) half-metallic SFMO when adding ASD. 


or an independence of ASD [14, 86] (see Fig. 6(a)). 

For all cases of electron correlations, Tc decreases 
with the inclusion of ASD (see Fig. 6(a)). The reduction 
is rather large for the GGA (AXc ~ 200 K) but becomes 
less in the case of half-metallic SFMO (AXc ~ 100 K). 
For the latter (Ujg = 2eV), the slope of reduction 
in Tc decreases for larger defect concentrations around 
10 at. %. This amount of ASD is a typical value found in 
real samples [85, 86], and the theoretical slope agrees 
rather well with the experimental observations or the 
previous Monte Carlo simulations [9]. In order to un¬ 
derstand the observed reduction, we have again a look 
at the magnetic coupling constants. They vary strongly 
as expected for the GGA or the Ujg = 2eV calculation 
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(see Fig. 7). We will mainly discuss the of half- 
metallic SFMO. By putting Fe at the Mo site (FeMo), a 
strong superexchange leads to a negative coupling con¬ 
stant (« —21meV) with the original Fe sites [10] while 
on the other hand the MoFe ions interact very weakly 
with the original Mo ions (not shown). The latter desta¬ 
bilizes the FiM ground state but at the same time the 
Fe-FeMo magnetic interaction should be able to counter¬ 
act these effects. Nevertheless, the FiM ground state 
collapses to an non-magnetic ordered ground state for 
more than 5at.% ASD and Tq drops down (Fig. 6(a)), 
since the AFM Mo-Fe becomes also reduced with increas¬ 
ing ASD (Fig. 7(b)). We want to note that for higher 
amounts of antisite disorder the superexchange interac¬ 
tions Fe-FeMo might stabilize a FiM order again with 
higher Tq or even dominate the magnetic coupling and 
form AFM order as observed in few experiments [88]. 

In contrast, the qualitative trend of Tq as a function 
of the oxygen content depends stronger on the electron 
correlation effects (see Fig. 6(b)). Tq becomes strongly 
reduced when using the GGA but it increases in half- 
metallic SFMO (Ujg = 2eV). For the latter, the mag¬ 
netic exchange coupling between Mo and Fe ions becomes 
again stronger and increases Tq by 15K/at.%. Electron 
correlation close to the fringe towards half-metallicity 
(U^g = 0.75 eV) showed a intermediate behavior for Tq 
- like for the total magnetic moment. All these obser¬ 
vations can be understood from the effect of the oxygen 
vacancies to the DOS and the magnetic exchange inter¬ 
actions. Both quantities vary with oxygen vacancies in a 
similar way as observed for increasing the Ujg parameter 
(cf. Figs. 2 and 5). This might be an artifact of the effec¬ 
tive medium description used for these calculations but it 
worked also very well for other oxide systems [23]. Thus, 
the oxygen deficiency leads to a localization of Fe 3d or¬ 
bitals and decreases by this the overlap integrals, which 
determine in a simple picture the exchange constants. We 
see for example that SFMO with = 0.75 eV becomes 
half-metallic when adding a certain amount of oxygen 
vacancies (minimum at ss 4at.% in Fig. 6(b)). 


D. The Curie temperature and electron doping 

Although we could find an estimation of the tendency 
of Tq as a function of the defect concentrations, the low 
absolute value might be attributed to various factors. We 
considered two possible mechanism, which might vary 
Tq as well, and could give another explanation for the 
discrepancy with the experimental Tq. 

On the one hand, electron doping might result from 
other defects, e.g., Navarro et al. [89] observed for SFMO 
an increasing Tq by substituting Sr with La. This was 
related with an increase of the density of states at the 
Fermi energy D(Ep). In theoretical calculations, such 
doping can be achieved by a shift of the Fermi energy 
Ep. For half-metallic SFMO, we indicate a variation by 
A Ep = 68meV in Fig. 8(a). The dashed line represents 


mm Fe e g Mo e g 
HI total DOS 5S53 Fe t 2g Mo t 2g 



E - E F (eV) E - Ef (eV) 

FIG. 8. LDOS of Fe and Mo around the Fermi energy calcu¬ 
lated with U Fe = 2eV for (a) a simulated electron doping and 
(b) with an exemplarily Umo = 2eV. The legend is provided 
above the plots. 


an n doping with additionally 0.15 electrons, which is 
related with an increase of D(Ep) and Tq (« 10 K). This 
results again from a stronger magnetic coupling between 
the Fe and Mo ions (see Fig. SI [71]). For a larger A Ep, 
stronger variations in the magnetic coupling constants 
beyond the next nearest neighbors (not shown) appeared 
whereby the Tq could not be increased further by electron 
doping. 

On the other hand, the electronic correlation effects 
were up to now only considered for the Fe d states. We 
found that the magnetic interaction and the Curie tem¬ 
perature are also very sensitive to the positions of Mo d 
states. The Mo d electrons of the majority spin channel 
are sharp and above the Fermi level, while in the minor¬ 
ity spin channel the d electrons are located at the Fermi 
level forming a relatively broad band (see Fig. 8(a)). We 
investigated exemplarily the influence of an additional 
Hubbard U parameter for Mo (Umo) hi a range of 0.5 eV 
to 4eV. A small shift of spin up electrons could sub¬ 
stantially change the Curie temperature. Since the DFT 
approach can not describe adequately unoccupied states, 
the positions of the d electrons of the majority spin chan¬ 
nel could be not correct. We did not observe any change 
of D{Ep) once we applied the Umo to all Mo d states. 
The main peak of the Mo dj g state retains its position 
above Ep but the unoccupied Mo states of the majority 
spin channel were shifted substantially to higher ener¬ 
gies. This led, e.g., to Tq = 250 K with Umo = 2eV (see 
Fig. 8(b)). 


VI. X-RAY ABSORPTION AND XMCD 
SPECTRA 


Finally, we want to return to the discussion of the va¬ 
lency of the Fe ions, which seems to be a mixed valency 
state following from experiments [28, 29, 33, 34] but is 
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FIG. 9. (Color online) The x-ray absorption spectra (Ref. [M3], 
open circles) at Fe L 2,3 edges as average of left and right 
circularly polarized light in SFMO measured at 10 K with 
5 T magnetic field compared with the theoretically calculated 
ones for Fe 3+ (dashed blue line) and Fe 2+ (solid red line). 


3+ in the ground state obtained by theoretical calcula¬ 
tions [10, 11, 72]. In a complex transition metal ionic 
compound such as the SFMO, the x-ray absorption and 
XMCD spectra at the L 2i 3 absorption edges can be used 
as fingerprints of the ground state. The experimentally 
measured Fe L 2i 3 absorption spectrum of the SFMO sin¬ 
gle crystal (Ref. [33], average of left and right circularly 
polarized light) displays at both absorption edges a weak 
lower-energy shoulder together with a doublet structure 
at the white line position with almost the same magni¬ 
tude (see Fig. 9). The complex fine structure of the Fe 
L 2j3 XAS is not compatible with a pure Fe 3+ valency 
state, since we observed only one peak - similar to [36] - 
instead of a double peak structure. In order to provide a 
quantitative description of the spectral features, we took 
into account both Fe valencies 3+ and 2+, separately, 
using a primitive unit cell. The theoretically calculated 
Fe L 2j 3 XAS agrees most closely with the experimental 
data by using those results in a linear combination of 
60% Fe 3+ and 40% Fe 2+ (see Fig. 9), which is opposite 
to the calculated proportion found in Ref. [33] using the 
framework of the ligand-field atomic formalism. 

Considering that the correct Fe valency in the SFMO 
ground state is 3+ as concluded in many theoretical stud¬ 
ies, the mixed valency state might follow again from lat¬ 
tice defects as speculated in [36]. We investigated the 
influence of the two types of defects in the tetragonal su¬ 
percell (see Fig. 1). For the antisite defects, we observed 
only the Fe 3+ solution. This might be connected with the 
high defect concentration modeled within the supercell. 
However, a single vacancy among twelve oxygen atoms 
has a more realistic concentration. Self-consistent calcu¬ 
lations in the tetragonal supercell produce the valency 
of the Fe ions being equal to 2.9+ and 2.4+ at the Fel 
and Fe2 sites, respectively. Therefore, the existence of 
the vacancy shifts the valency of the nearest Fe ion (Fe2) 



FIG. 10. (Color online) The x-ray absorption spectra 
(Ref. [33], open circles) at Fe L 2 ,3 edges in SFMO measured 
with left (<r + , top panel) and right circularly polarized light 
(<r - , middle panel) measured at 10 K with 5T magnetic field 
and XMCD experimental spectrum (Ref. [33], lower panel) 
compared with the theoretically calculated spectra with an 
Fe ion far away (dashed blue line) and near to an oxygen 
vacancy (solid red line). 


towards 2+. This valency change could be also observed 
in the CPA calculation of the DOS including oxygen va¬ 
cancies (see Fig. 2). 

Indeed, the full explanation of the experimental spec¬ 
tra is only possible by taking into account these crystal 
imperfections. The Fe L 3 x-ray absorption spectrum for 
left circularly polarized light (cr + ) possesses four major 
fine structures a, b, c and d (see Fig. 10). We found 
that the calculations for the ideal crystal structure with 
the Fe 3+ ground state solution (dashed blue line) pro¬ 
vides the x-ray absorption intensity er + only at the major 
peaks 6 , c, and d and do not reproduce the low energy 
shoulder (peak a) as well as the low energy peak e at 
the L 2 edge. The calculations with the Fe 3+ solution 
produce only one high energy peak structure in the L 3 
<r~ spectrum (middle panel). However, the experimen¬ 
tal measurements exhibit a double-peak structure. The 
x-ray absorption from the Fe2 atoms with the oxygen va¬ 
cancy (solid red line) contributes to the low energy peak 
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FIG. 11. (Color online) The experimental x-ray absorption 
spectra (Ref. [ 35 ], open circles) at O K edge (top panel) in 
SFMO and experimental XMCD spectra (lower panel) mea¬ 
sured at T = 20 K with B = 1.1T compared with the theo¬ 
retical simulations carried out for the ideal crystal structure 
(solid red line) and with an oxygen vacancy (dashed blue line). 



Energy (eV) Energy (eV) 

FIG. 12. (Color online) The experimental x-ray absorption 
spectra (Ref. [33], open circles) at Mo L 2 ,3 edges (top panels) 
in SFMO measured at 10 K with 5T magnetic field and ex¬ 
perimental XMCD spectra (lower panels) compared with the 
theoretical simulations. 



Energy (eV) 

FIG. 13. (Color online) The experimental x-ray absorption 
spectra (Ref. [33], open circles) at Mo M 2,3 edges [90] (top 
panels) in SFMO and experimental XMCD spectra (lower 
panels) compared with the theoretical simulations. 

of a~ absorption (see Fig. 10). The relative intensity of 
the peaks depends on the relative concentration of the 
Fel and Fe2 ions in SFMO, in other words, the concen¬ 
tration of defects such as oxygen vacancies. It is similar 
for the oxygen K edge. The calculations including an 
oxygen vacancy are in better agreement with the exper¬ 
imental measurements in the x-ray absorption as well as 
in the XMCD (see Fig. 11). 

In contrast, the XAS and XMCD spectra at the Mo 
L 2 ,3 and M 2 ,3 edges are less sensitive for the crystal de¬ 
fects. For both edges the agreement with the experimen¬ 
tal measurements is quite good (see Figs. 12 and 13) [ 3]. 


VII. CONCLUSIONS 

We performed an extensive study of the electronic and 
magnetic properties of Sr 2 FeMoOg with the main focus 
on the influence of defects. In order to agree with ex¬ 
periments and to describe the half-metallicity of SFMO, 
we had to apply an U e g parameter on the d states of Fe. 
The band gap in the upper spin channel opened already 
at a small value of Ujg = 1 eV. Otherwise, the system 
remained metallic. The best agreement with available 
experimental data for XAS, XMCD, To and DOS was 
achieved at Ujg = 2 eV. As already discussed in previous 
works, the main contribution at the Fermi energy orig¬ 
inates from the hybridized t 2g states of the Mo and Fe 
ions. 

Only the consideration of a magnetic moment at the 
Mo site stabilized the ferrimagnetic ground state and 
yielded a reasonable To Both kinds of lattice defects 
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ASD and V va c — reduced the total magnetic moment, 
which is in good agreement with the experimental stud¬ 
ies. On the other hand, we found several mechanisms, 
which affected Tq. If the measured sample was doped 
with additional electrons or electron correlation effects 
were taken into account also for the Mo ions, Tq could 
be increased. Antisite disorder might have had the op¬ 
posite effect and reduces Tq. The most important result 
is the prediction of an increasing Tq when considering 
oxygen vacancies as well. These vacancies seem to be 
also responsible for the specific magnetic dichroism spec¬ 
tra for the Fe 1/2,3 edges observed in SFMO. We obtained 
only good agreement with the measured spectra including 
oxygen deficiency. It might be possible to estimate the 
amount of oxygen vacancies via fitting the experimental 


XAS/XMCD spectra with a particular mixture of theo¬ 
retical spectra for defect-free and oxygen deficient SFMO. 
Finally, we can can conclude that varying defect concen¬ 
trations could be in principle one reason for the observed 
differences in the magnetic properties of thin films or bulk 
samples discussed in the introduction. The appearance 
of defects might vary Tq in the order of ~ 15K/at. %. 
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Supplemental Material 

e. Varying magnetic coupling constants by electron doping In order to model electronic doping on a 
simple level, we varied for the calculation of the magnetic exchange parameters Jij the Fermi energy, Ep, by several 
meV without recalculating the self-consistency loop for the density of states (DOS). The corresponding nearest and 
next nearest magnetic exchange parameters are plotted in Fig. SI. The antiferromagnetic (AFM) coupling between 
the Fe and Mo ions became again stronger. This increased also the critical temperature, Xc, by few K. 



FIG. SI. Nearest and next nearest Jo) e " A for the coupling between Fe-Fe and Fe-Mo as a function of non-selfconsistent variation 
of Ep. X is either Fe or Mo and i indicated the shell number. 

f. Electronic structure and electronic correlations in SFMO In order to determine an appropriate Ufg 
value, we aimed to match the band gap of the majority spin channel with experimental data. The most relevant 
variations close to the Fermi energy happened for the Fe and Mo d states (see DOS obtained from GGA calculations 
in Figs. S2(a) and S3(a)). The measurement of the photoelectron spectrum (PES) by Saitoh et al. (SlJ allowed an 
estimation of the position of the Fe e g states (« 1.3 eV) assuming their comparison with a theoretical DOS calculation 
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FIG. S2. Atomic and orbital resolved DOS for the Fe d states: (a) GGA calculation and (b) GGA+I7. In (b) the DOS is 
color-coded as a contour plot. The dashed line in the left hand side panel indicates the peak position in the experimental PES 
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FIG. S3. Atomic and orbital resolved DOS for the Mo d states: (a) GGA calculation and (b) GGA+LL In (b) the DOS is 
color-coded as a contour plot. 

is valid. When increasing the Ujg parameter, the main changes at the Fermi energy appeared for the Fe e g states (see 
Fig. S2(b), left hand side). They were pushed towards lower energies while all other states in the spin up channel did 
not cross Ep (see Figs. S2(b) and S3(b)). Thus, a band gap in the majority spin channel opened up at <1 leV. 
In contrast, hybridizing states from the Fe and Mo dj g states remain in the minority spin channel at Ep, which is why 
we observed half-metallic SFMO. 


[SI] T. Saitoh, M. Nakatake, A. Kakizaki. H. Nakaiima, O. Morimoto, S. Xu, Y. Moritomo, N. Hamada, and Y. Aiura, Phys. 
Rev. B 66, 035112 (2002). 























